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We derive the contributions of spin-orbit and spin-spin coupling to the gravita- 
tional radiation from coalescing binary systems of spinning compact objects. We 
calculate spin effects in the symmetric, trace-free radiative multipoles that deter- 
mine the gravitational waveform, and the rate of energy loss. Assuming a balance 
between energy radiated and orbital energy lost, we determine the spin effects in the 
evolution of the orbital frequency and orbital radius. Assuming that a laser interfer- 
ometric gravitational observatory can track the gravitational-wave frequency (twice 
the orbital frequency) as it sweeps through its sensitive bandwidth between about 10 
Hz and one kHz, we estimate the accuracy with which the spins of the component 
bodies can be determined from the gravitational-wave signal. 
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I. INTRODUCTION 



The ability of laser interferometric gravitational-wave (GW) detectors such as the US 
LIGO to extract useful astrophysical information about coalescing binary systems of com- 
pact objects from the observed gravitational-wave signals depends upon highly accurate 
theoretical models for the decaying orbital evolution. Detection and study of the character- 
istic "chirp" waveform emitted by such systems involves a matched filtering technique using 
a theoretical template that is a function of parameters of the source [|T|. As has been pointed 
out by Cutler et al. [Q, matching the theoretical variation of the GW frequency with the 
data as the frequency sweeps through the detector bandwidth between 10 and 1000 Hz may 
be the most promising way to obtain estimates of the masses and spins of the bodies. To this 
end, one needs a theoretical formula for the evolution of the GW frequency that includes 
all dependence on the two masses and spins, and that is sufficiently accurate that errors in 
the model lead to cumulative phase errors no larger than 2ti radians over the potentially 
thousands of cycles observed in the frequency bandwidth of interest (larger phase errors 
substantially reduce the signal-to-noise of the matched filter). 

Cutler et al. p[ have suggested that a model for the frequency evolution that is based on 
a post-Newtonian approximation may be inadequate unless the approximation is carried to 
extremely high orders beyond the lowest-order Newtonian and post-Newtonian corrections 
that have been derived to date PJ^. In addition, they point out that the effects of spin-orbit 
and spin-spin coupling of the bodies have not been derived, even in the first post-Newtonian 
approximation. It is the purpose of this paper to remedy the latter situation. 

We derive the spin-orbit and spin-spin contributions to the symmetric trace-free (STF) 
multipoles that enter the gravitational waveform for two-body systems, using a multipole 
formalism developed by Blanchet, Damour and Iyer (BDI) [^^. In addition to explicit spin 
terms, the multipole formalism involves several time-derivatives of multipole moments that 
lead to expressions involving the two-body acceleration; to evaluate these moments consis- 
tently, we use post-Newtonian equations of motion that also include spin terms. Using these 
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multipoles, we then evaluate the rate of energy loss from the binary system. Specializing to 
a circular orbit, we calculate the rate of decrease of orbital radius and the rate of change of 
orbital frequency, assuming energy balance and using an expression for energy that includes 
spin terms. Finally, we estimate the size of spin-orbit and spin-spin effects for coalescing 
binary systems of neutron stars and/or black holes on the accumulation of GW phase in a 
LIGO-type detector. 



II. EQUATIONS OF MOTION 

Equations of motion for bodies of arbitrary mass and spin have been developed by numer- 
ous authors (for reviews and references see For our purposes, only the Newtonian, 
spin-orbit (SO) and spin-spin (SS) terms will be needed; post- and post-post-Newtonian 
non-spin terms have been considered elsewhere 0. Although the inclusion of spin implies 
that we are treating extended bodies, we ignore tidal and quadrupole-coupling effects (even 
though rotationally induced quadrupole effects are proportional to (spin)^); for binary sys- 
tems containing neutron stars or black holes, these are expected to be small until the very 
latest stage of inspiral and coalescence, except possibly for rapidly rotating Kerr black holes 

By eliminating the center of mass of the system, we convert the two-body equations of 
motion to an effective one-body equation of motion given by (G = c = 1) 

TTl 

a = ai - aa = — + aso + ass, (1) 



where 



aso = ^ {6n [(n x v) ■ (C + 0] " [vx (4C + 30] + 3r [fix (2C + 0]} , (2a) 

ass = {n (Si ■ S2) + Si (fi ■ S2) + S2 (fi ■ Si) - 5n (n ■ Si) (fi ■ S2)} , (2b) 

jj,r^ 

where m = mi + m2, /i = mim2/m, r = |x|, x = xi — X2, n = x/r, v = Vi— V2, C = S1-I-S2, 
^ = ("^2/'"^i)Si + (mi/m2)S2, and r = (fi ■ v). 



It is useful to note that spin-orbit and spin-spin effects are of order {R/r)vv and {R/r)'^v'^, 
respectively, compared to the Newtonian acceleration, where R and v denote the charac- 
teristic radius and rotation velocity of each body; these terms thus are formally of first- 
post-Newtonian order. However, for compact bodies, R is of order of m or a few times m, 
and V could be of order unity, so that in practice these terms can be considered to be of 
post^/^-Newtonian and post-post-Newtonian order, respectively (indeed they are so denoted 
in i). 

Equations for the precession of the spins caused by spin-orbit (geodetic precession) and 
spin-spin (Lense-Thirring) couplings can also be written down in terms of relative 

coordinates. The relevant equation for our purposes is 

C = ^LnX + 1^)+^ [(A • S2) (n X Si) + (n ■ Si) (fi x S2)] . (3) 

where L^r = /ux x v denotes the Newtonian orbital angular momentum. Note that the 
precession of the spins is one order higher than the spins themselves. Therefore we may 
neglect this precession when evolving the equations to lowest order; Eq. (^ is only needed 
to verify that total angular momentum is conserved for the system, apart from that radiated 
away [p]) . 

These equations of motion can be derived from a generalized Lagrangian which is a 
function of the relative position, velocity, and acceleration, given by 

1 TTl 

C = -fiv'^ + fi h Cso + Css, (4) 

2 r 

where 

Cso = I- [v- (a X 0] + 24 {v [XX (C + 0]} , (5a) 

Css = ^ {(Si • S2) - 3 (n • Si) (fi ■ S2)} . (5b) 

Here the Euler-Lagrange equations are dC/dx"^ — dpi/dt = 0, where = dC/dv^ — s*, and 
= dC/da^. It is understood that wherever the acceleration appears in higher-order terms 
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in the Euler-Lagrange equation, one substitutes the lower-order equation of motion (see |T2 
for discussion of acceleration-dependent Lagrangians). 

The relative Lagrangian is invariant with respect to time translations so that there exists 
a conserved energy, given hj E = p*f * + s*a* — C Evaluating this expression we obtain 

X ITl 

E = -fiv^ - /i- + Eso + Ess, (6) 
2 r 

where 

i?5o = ^{(fix v)-a, (7a) 

Ess = {3 (fi ■ Si) (n • S2) - (Si • S2)} . (7b) 
We can define the total angular momentum as L = + (x x p) + (v x s), with the result 

L = C + + Lso, (8) 

where 

^{-[hx (fix [2C + m - I [vx (v X e)]| , (9) 
m I r 2 J 



■'SO 



where there is no spin-spin contribution to L. Using the equations of motion and spin 
precession, Eqs. (P and (^), it is straightforward to show explicitly that, to post-Newtonian 
order, E = h = 0. 

III. EVALUATION OF BDI MULTIPOLES 

The radiative energy loss of the system can be expressed in terms of symmetric and 
trace-free (STF) radiative multipole moments (see |jl3|] for a review). For the accuracy we 
require the energy loss rate is given by 

dE _ _lf(p (3) ^ (4) (4) 16 (3) (3) ^ 
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where lij and lijk are the STF "mass" quadrupole and octopole moments, and Jij is the 
STF "current" quadrupole moment, and (n) over each moment denotes time derivatives. 

We evaluate the radiative multipoles using the formalism developed by Blanchet, Damour 
and Iyer [^,|| (BDI). We restrict ourselves to the case of two well-separated, approximately 
spherically symmetric, rotating compact objects whose structure is given by that of a per- 
fect fluid. Following our previous post-Newtonian approach []3|,[l4|, we choose the following 
definition for the center of mass of each body 



x\ = — j x'p* (x) 

A 



1 + + 



d^x, (11a) 



mA = J p ix 

A 



1 + + 



£x, (lib) 



where p* = p(l + |f ^ + 3U) is the so-called "conserved density" [|15|, with p the local mass 



density, v the velocity, and U the Newtonian gravitational potential; 11 is the specific internal 
energy density, v\ = v'' — v\, v\ = dx^/dt, and f/4 is the Newtonian potential produced by 
the A-th body itself. 

However, when spin effects are to be included, there is a subtle difference between the 
center of mass defined above, and the center of mass used in the equations of moton (p. 
The latter defines the center of mass world line x^'^ of each body using a so-called "spin 
supplementary condition" (SSC), given by Sa'^'^uai, = 0, where Ua is the four- velocity of the 
center-of-mass world line, and 

Sa"" = 2 J {X^" - XA^^')T''^''d^X, (12) 

where r'^^ denotes the stress-energy tensor of matter plus gravitational fields, satisfying 
t'^" I, = 0, and square brackets around indices denote antisymmetrization. Note that the 
spin S of each body is defined by S"^ = ^eijkSA- 

It is then simple to show that, evaluating the BDI multipoles using the center of mass 
definition Eq. (|lla|) , and then making the transformation 
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X A + 



2mA 



(VA X Sa)' 



(13) 



we can convert all expressions to the center of mass defined by the SSC of the equations of 
motion. Because the correction is of post-Newtonian order, it needs to be made only in the 
Newtonian-order multipoles. 

A useful check of the consistency of our approach comes from evaluating the mass dipole 
moment of the system, given by Eq. (Al6b) of Ref. P], and making the transformation 
(O). The result is /' = Yj\f^Ax\ + (va x Sa)*)-- The original two-body equations of 
motion from which Eq. ([l|) is derived |P] then imply that P = 0, as expected (uniform 
motion of the system's center of mass). 

Evaluating the other multipoles, transforming them to the SSC of our equations of motion 
using Eq. (|13]), and then transforming them into relative coordinates, we obtain 



rii ( i 8 
1 = jji\x X-' j + -r] 

jijk 



X* (v X 



STF 



v' (x X 



STF 



STF 



( 6m\ 


X* (x X v)-' 


STF 3 




+ 2 


\ m ) 







, STF 



(14a) 
(14b) 
(14c) 



where 5m = mi — m2, rj = fi/m, and a = (m2/m)Si — (mi/m)S2. Note that, to post- 
Newtonian order, there are no explicit spin-spin contributions to the multipoles. 



IV. ENERGY LOSS AND INSPIRAL OF CIRCULAR ORBITS 



Taking time derivatives of Eqs. (14), substituting the equations of motion where appro- 
priate, and substituting the results into Eq. ([10|) , we obtain the energy loss rate 



dE 8 mV^ 



dt 



+ 



15 

fi X V 



mr 
1 



C 27r' - 37v^ - 12- 



m 



^ 51r^ - 43f ^ + 4- 



, m 



3(Si ■ Sa) (47i;2 - SSr^) - 3(n • Sin ■ 83) (l68v^ - 269f^'^ 



2mfir'^ 

+71(v ■ Siv ■ S2) - 171r (v Sifi ■ S2 + n ■ Siv ■ 83)]} . (15) 



We now restrict ourselves to the case of nearly circular orbits, by which we mean orbits 
whose inspiral timescale due to gravitational radiation energy loss is long compared to an 
orbital timescale, and in which r = constant+ small, periodic perturbations due to spin 
couplings (when the spins are orthogonal to the orbital plane, the spin perturbations are 
constant). To simplify the discussion, we then take an angular average of all quantities over 
an orbit. We then have f = r = n- v = 0, = r^f2^, where Q is the angular frequency, 
Ln = fir'^Qh, where L is a unit vector orthogonal to the orbital plane. The equation 
n ■ a = f — rfi^ yields the following orbit-averaged relationship for a circular orbit 



These conditions lead to the following results for the energy and energy loss rate for a 
circular orbit, averaged over an orbit, expressed in terms of Q: 




2 fim^/'^ 



(Si • S2 - 3L ■ SiL • S2) 



(16) 




(Si ■ S2 - 3L ■ SiL -82) , (17a) 




(17b) 



The result for the evolution of orbital frequency is 




(18) 



The energy and loss rate can also be expressed exclusively in terms of r using Eq. (jTB]), 



and an equation for the rate of inspiral obtained. The result is 




(19) 
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V. DISCUSSION OF RESULTS 



A signal template whose frequency evolution is given by Eq. (|T^ can be said to match 
the signal if the accumulated phase in GW from the time the signal enters the detectors' 
sensitive bandwidth to the time it leaves it differs from that of the template by less than 
2-71 radians. Thus one can obtain a crude estimate of the accuracy with which a parameter 
characterizing the template can be determined by finding that change in the parameter 
that leads to a change of 27r in the accumulated phase Such estimates are only crude 
and probably optimistic, because they do not take into account signal-to-noise issues or 
correlations among multiple parameters. 

The accumulated phase in gravitational waves is given by 

C f P f 

^Gw = '2n fdt = 2 {QyQ)dQ/Q, (20) 
where / = fi/vr is the GW signal frequency. Integrating the SO and SS terms in Eq. (|18|), 



and ignoring precession of the spins or of the orbital angular momentum, we can estimate 
the accuracy of determination of the spin parameters. For example, for two equal-mass 
neutron stars, we find, for the spin-orbit terms, 

5/3 / f \ 2/3 



^ ^ 38 ( V (J^\ (Jjll.) 



(21) 



where rj^s and P^s are the neutron-star radius and rotation period (suitably averaged over 
the two stars), and /in is the GW frequency entering the detector bandwidth. For a 10 Mq 
black hole with spin = rriiai, where ai < mi is the Kerr parameter, and a 1.4 Mq neutron 
star, the quantity ai/mi be determined to an error A(ai/mi) ~ 0.016 for ai/mi > 0.03, 
while for ai/mi <^ 0.03, the spin of the companion neutron star cannot be found to better 
than a factor of four, for a 10 ms rotation period. For two 10 Mq black holes, the error 
in the net "Kerr parameter per mass" of the system projected orthogonal to the orbital 
plane, A(ai/mi + a2/m2)z, can be estimated to be about ±0.1. However, Monte-Carlo 
studies indicate that strong correlations between post-Newtonian terms (dependent on /i) 
and spin-orbit terms are likely to weaken these estimates substantially p|. 



Spin-spin terms, on the other hand, have neghgible effect on the accumulated phase for 
most systems of interest. Only for two extreme Kerr black holes are these terms discernable; 
for two IOMq black holes with both spins orthogonal to the orbital plane, we estimate 
A{SiS2)/\SiS2\-0A. 
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